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O ! ABSTRACT 

O 

^ ■ The Gil. 11-0. 12 infrared-dark cloud has a filamentary appearance, both in 

extinction against the diffuse infrared emission of the Galactic plane and in emis- 
sion at 850/ini. We use a novel computational technique based on an advanced 
genetic algorithm to explore thoroughly 3 different models of self-gravitating, 



pressure truncated filaments and to constrain their parameters. Specifically, the 
models tested are the non- magnetic Ostriker (1964) model, a generalized ver- 
sion of the magnetic Stodolkiewicz (1963) model, and the magnetic Fiege & 
Pudritz (2000a) model. Previous results showed that Gil. 11-0. 12 has a much 
steeper ~ r~ 4 radial density profile than other filaments, where the density varies 
q ■ approximately as r -2 , and that this steep density profile is consistent with the 

Ostriker (1964) model. We present a more complete analysis that shows that 
the radial structure of Gil. 11-0.12 is consistent with regimes of each of these 
models. All of the magnetic models that agree with the data are threaded by a 
dominant poloidal magnetic field, and most have dynamically significant fields. 
Thus, Gil. 11-0. 12 is an excellent candidate for radial support by a magnetic field 
that is predominantly poloidal. We predict the polarization patterns expected for 
both magnetic models and show that the two magnetic models produce different 
polarization patterns that should be distingished by observations. 

Subject headings: ISM: clouds — magnetic fields — MHD 



1. Introduction 

Infrared-dark clouds (IRDCs) are a class of molecular clouds seen in extinction against 
the diffuse infrared emission of the Galactic plane. About 2000 of these objects were discov- 
ered by the Midcourse Space Experiment (MSX) satellite in spectral bands from 7 — 25/im 
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(Egan et al. 1998). They reside at distances of 2 — 8 kpc typically (Egan et al. 1998; Carey 
et al. 1998), and most are located roughly toward the Galactic centre. The IRDCs have 
unusual physical characteristics compared to nearby examples of molecular clouds. Egan 
et al. (1998) found that IRDCs appear to have properties intermediate between massive 
molecular cloud cores (n ~ 10 4 — 10 6 cm" 3 ; T » 30 — 100 K) and isolated dark clouds 
(n ~ 10 2 — 10 4 cm~ 3 ; T 8 — 15 K), having densities typical of molecular cloud cores, but 
temperatures typical of isolated dark clouds. Many IRDC filaments appear to be at a rela- 
tively early phase of evolution, with fewer signposts of star formation than nearby molecular 
clouds. 

The IRDCs and nearby molecular clouds share the common feature that a large fraction 
of both populations are filamentary, but the IRDCs are often less fragmented. Relatively 
"pristine" IRDC filaments have simpler density structures than fragmented, star-forming fil- 
aments, which makes them good candidates for detailed comparison with theoretical models. 
One drawback is that IRDCs are more distant than many other filaments and are therefore 
less well resolved. 

The Gil. 11-0. 12 cloud studied in this paper is a good example of a filamentary IRDC, 
which was observed by SCUBA in 850/xm emission Carey et al. (2000), following its discov- 
ery by the MSX survey. Johnstone et al. (2003) (hereafter Paper I) fit a non-magnetic 
equilibrium filament model to the radial structure of two regions of Gil. 11-0.12. This pa- 
per presents a more thorough analysis, which compares the SCUBA 850/zm emission with 
3 existing models of filamentary clouds in detail. These models are the pressure-truncated 
non-magnetic Ostriker (Ostriker 1964) model, a generalized Stodolkiewicz (1963) (hereafter 
GS) helical field model, and the Fiege & Pudritz (hereafter FP) helical field model. We 
discuss these models in the next section. Note that the Ostriker model is the zero magnetic 
field limit of both the GS and FP models, but this limit is sufficiently important to warrant 
separate treatment. 

This paper uses a novel genetic algorithm-based technique to explore the parameter 
space of each model completely and to fit solutions to the data. This is the first astronomical 
research application for a new multi-object genetic algorithm ("Ferret") developed by one of 
the authors (J.F.). The technique is discussed in Section 5. 

2. Theoretical Models of Filamentary Clouds 

The goal of this paper is a thorough comparison of the Gil. 11-0. 12 IRDC with 3 differ- 
ent theoretical models of filamentary molecular clouds. These models have different magnetic 
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field configurations, or no magnetic field at all in the case of the Ostriker model. However, 
they share the common features of self-gravity, external pressure truncation and isother- 
mality, meaning that the velocity dispersion is constant everywhere within the cloud. This 
last feature is somewhat restrictive because the velocity dispersion is well-known to increase 
with size scale in molecular clouds, and the equation of state does affect the radial density 
structure. However, the effective equation of state of the gas is not known, and we restrict 
our analysis to existing, classical models of filaments. 



2.1. Non-Magnetic Filaments: The Ostriker Model 

The simplest equilibrium model of a filamentary cloud is that of a self-gravitating, non- 
magnetic, isothermal filament. Ostriker (1964) studied such clouds, and found an exact 
solution for the density structure 1 

P= 2, (1) 

[l + rV(8r 2 )] 2 ' 

where r is the core radius of the filament. The core radius is defined by 

where a is the 1-dimensional velocity dispersion of the gas, assumed to be constant and pos- 
sibly including both thermal and non-thermal contributions, G is the gravitational constant, 
and p c is the central density along the axis of the filament. Note that the density varies as 
~ r~ 4 at large radii r ^> r . 

In the absence of external pressure, the Ostriker solution extends to infinite radius, and 
the mass per unit length, which we refer to as the line mass for brevity, is exactly 

m vir = q ■ (3) 

This quantity is denoted m vir , since it also represents the virial line mass in the virial equation 
4 below. 

An Ostriker filament is truncated at finite radius if it is bounded by an external pressure. 
In this case, the filament obeys the virial equation 

— = 1 " — , (4) 

TYl v i r (P ) 



1 This solution is equivalent to the non-magnetic limit of the magnetic model discovered by Stodolkiewicz 
(1963), discussed in §2.2. 
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where m is the line mass, Ps is the external pressure, and (P) is the mean (volume weighted) 
internal pressure defined as follows: 

(P) = £™ (5, 

(Fiege & Pudritz 2000a) (hereafter FP1). Generally, we use calligraphic characters to refer 
to quantities defined per unit length. Thus, V = rtRs 2 is the cross-sectional area of the 
filament (volume per unit length), in which Rs is the radius of the filament where it is 
truncated by the pressure of the external medium. 

Thus, the line mass of the Ostriker filament is reduced if the filament is bounded by a 
finite external pressure. Filaments with a line mass exceeding m vir cannot be represented 
by an hydrostatic equilibrium model and require magnetic support against self gravity. 

We parameterize the pressure truncation radius Rs of the Ostriker model by the con- 
centration parameter 

C - log 10 (^) • (6) 

Thus, the Ostriker model is characterized by a single dimensionless parameter C, plus two 
parameters p c and a, which determine the density scale and core of the solution via equation 
2. 



2.2. The Generalized Stodolkiewicz Model 

Stodolkiewicz (1963) discovered an analytic solution for isothermal self-gravitating fila- 
mentary clouds threaded by a poloidal (B z ) magnetic field. Stodolkiewicz assumed that the 
magnetic pressure is proportional to the gas pressure, so that 

P = p = — 52~ = C0USt ( ? ) 

within the cloud. 

Magnetostatic models always have the freedom to assume how the mass is distributed 
relative to the poloidal magnetic flux surfaces. The two most common assumptions are that 
j3 is constant, or that the mass (per unit length for a filament) and magnetic flux enclosed 
by each poloidal magnetic flux surface are proportional. The Stodolkiewicz model is of the 
former type, while the FP model, discussed in §2.3 below, is of the latter. 

The Stodolkiewicz model was generalized by Nakamura, Hanawa and Nakano (Naka- 
mura et al. 1993) to include a toroidal (B^) magnetic field component, so that the field lines 
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are helical in general. The radial density structure is identical to that of the Ostriker model 
(equation 1), except that the core radius r is replaced by r 0tS tod, which is modified by the 
magnetic field: 

1 + cos 2 6T 



r 0,Stod 



1 + 



2(3 



(8) 



where is the pitch angle of the magnetic field, relative to the axis, at radial infinity. The 
poloidal and toroidal components of the magnetic field are given respectively as follows: 



B z = B z (0) 



-3/2 

co ^[^I^l 11 + 33—1 ( 9 ) 




"o = B,(0) ^— sine[U - ] • (10) 

where B z (0) is the magnetic field at the radial centre of the filament, given by 



b.(o) = J*p. (ID 

The pressure-truncated GS model is characterized by 3 dimensionless parameters C, (3 
and 9, where C is the same concentration parameter defined by equation 6, plus plus two 
parameters p c and a, which determine the density scale and core radius exactly as in the 
Ostriker solution. Note that there exists an interesting and exact degeneracy in the density 
structure of the Stodolkiewicz (1963) model. All solutions characterized by 

(1 + cos 2 9)/2p = k (12) 

form a 1-parameter set having exactly the same core radius and density structure, controlled 
by the constant k, although they may have different magnetic field strength and geometry, 
which are controlled by 9 and (3. This paper fits models to the surface density structure of 
Gil. 11-0. 12 only. Thus, we consider only k in our analysis. 

The density structure of the GS model tends to that of the Ostriker model in the non- 
magnetic limit (3 — > oo. Note that r^stod > fo for any choice of (3 and 9, which indicates 
that the magnetic field always offers net magnetic support. Thus, the GS helical field model 
is always dominated by the poloidal magnetic field. This is true even for the special case 
9 = 7r/2, where the magnetic field tends to a purely toroidal field as r — > oo. Even in this 
limit, the poloidal field dominates globally over the toroidal field and provides net support. 
The GS model differs in this regard from the FP model, which may be dominated by either 
the poloidal or toroidal field. 
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2.3. The Fiege & Pudritz Model 



The primary difference between the GS model and the FP model (Fiege & Pudritz 
2000a) is that they make different assumptions regarding the magnetic flux loading. The 
FP model assumes uniform flux to mass loading, whereas the GS model assumes that the 
plasma f3 parameter is constant (see (§2.2). The assumption of a constant magnetic flux to 
mass ratio in the FP model implies that that all magnetic field lines are loaded with equal 
mass. Thus all field lines are equivalent, and there is no possibility of exciting magnetic 
interchange instabilities. Also, any filament with this property could be assembled from a 
cloud of uniform density and magnetic field by radially compressing and twisting flux tubes 
under the condition of magnetic flux freezing. 

The magnetic flux to mass ratios T z and are defined as follows: 



which determine the flux to mass loading of the poloidal and toroidal field components 
respectively. We refer the reader to FP1 for an illustration of T z and T^. FP1 found that 
the virial properties of the filaments that they studied could be described by the isothermal 
FP model with 



In contrast to the Ostriker and GS solutions, the density and magnetic field of the FP 
model cannot be written down in closed form. One must integrate numerically the equations 
of self-gravitating magnetostatic equilibrium, which are given in FP1. However, the magnetic 
generalization of equation 4 is only slightly more complex: 



where Pg, (P), m, and m vir are defined as above, and the remaining terms are defined as 
follows. The gravitational binding energy per unit length is given by 




(13) 



5<r^ < 25 
r, < 8. 



(14) 




(15) 



>V = - m 2 G, 



(16) 



and the total magnetic energy per unit length is given by 




(17) 
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In the equations above, subscript S refers to quantities evaluated at the cylindrical surface 
of the cloud. 

The assumption of constant magnetic flux to mass ratio makes the FP model somewhat 
more general than GS model, in the sense that the magnetic energy term M. in equation 
15 may be either positive, negative or zero, corresponding to net support, net magnetic 
confinement or magnetically neutral solutions, respectively. The magnetic energy term in 
the GS model is always positive, or zero in the special case where the solution tends to the 
non-magnetic Ostriker limit discussed in §2.1. 

The sign of M. in equation 15 determines whether the magnetic field provides net mag- 
netic support against self-gravity [M. > 0) or confinement (M. < 0). The former case 
(M. > 0) is the result of a dominant poloidal magnetic field, for which the field lines provide 
a radially expansive pressure. The latter case (Ai < 0) results from a dominant toroidal 
magnetic field, which provides a net binding energy. The FP model applies equally well to 
either of these regimes. FP1 tabulated virial data that allowed them to estimate Ps/{P) 
and m/m vir for a number of filamentary clouds, which provided an estimate the magnetic 
energy term A^/|W|. They found that A^/|W| < for most of their sample, which is only 
possible for filaments dominated by the toroidal field. 

Solutions with flux to mass parameters in the ranges given by equation 14 may be 
dominated by either the poloidal or the toroidal field component. Toroidal- field dominated 
filaments are characterized by density profiles that are more centrally peaked than the Os- 
triker model, but fall off more gradually with radius. FP1 found numerically that the density 
falls off as r~ L9±ai at large radius, which is in agreement with the observed density profiles 
of several filamentary clouds (Alves et al. 1998, Lada et al 1999, Johnstone & Bally 1999). 
Fiege (2003) showed analytically that this shallow density profile is due to the dominance 
of the toroidal field at large radii (r ^> r ). Steeper density profiles are possible in the 
regime dominated by the poloidal field in the FP model. FP1 studied this regime, but they 
focused on the toroidal-field dominated regime of the model, which was more consistent with 
observational data compiled from the literature. 

Poloidal field and toroidal-field dominated regimes of the FP model are separated by a 
magnetically neutral regime of the model, defined by M. = 0, where the expansive pressure 
of the poloidal field is exactly balanced by the radial confinement due to the toroidal field. 
The virial equation 15 cannot distinguish solutions in this neutral set from the non-magnetic 
truncated Ostriker solution. This illustrates a fundamental degeneracy of the FP model. 
It is generally true that any choice of A^/|W| and surface pressure corresponds to a one 
parameter set of models, where T z and are critically balanced to yield the same A^/|W|. 
This degeneracy occurs because B z and enter the magnetic energy term of the virial 
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equation (equation 17) with opposite sign, so that poloidal and toroidal terms can balance 
each other. In general, FP models with indistinguishable virial properties can still have 
observable differences in their internal density structure. This study is based on detailed 
models of the internal density structure, which might distinguish between such differences. 

To summarize, the FP model requires 3 dimensionless parameters, namely the concen- 
tration parameter C and two magnetic flux to mass ratios T z and T^, which correspond 
respectively the flux to mass ratios of the poloidal and toroidal magnetic field components. 
These dimensionless parameters determine the shape of the density, pressure, and magnetic 
profiles. Two dimensional parameters p c and a are also required, which determine the den- 
sity scale and core radius as in the other models. C, p c and a are as defined for the Ostriker 
and GS models. Both the GS and FP models assume that the field may be helical in general, 
although a pure poloidal field is allowed by both models. 

3. Analysis of the Gil. 11-0. 12 Filament 

Figure 1 shows an 850/im map of the Gil. 11-0.12 filamentary infrared-dark cloud. The 
observational details and data reduction are described in Paper I. At an estimated distance 
of 3.6 kpc (Carey et al. 1998), the Gil. 11-0. 12 filament is ~ 22 pc long and has a width 
of ~ 1 pc. The filament is curved and fragmented into several cores, but none of these 
cores appears to be forming stars. Thus, the fragmentation and subsequent evolution is 
probably in a less advanced stage than in other well-studied filaments such as the Orion 
"integral-shaped" filament (Dutrey et al. 1993; Johnstone & Bally 1999). 

The first step of our analysis was to locate accurately the central "spine" of the filament. 
We began with an initial guess by graphically locating 10-20 points along the ridge covering 
almost the entire length of the filament. We smoothed this curve slightly using a low-order 
windowed polynomial fit, and constructed orthogonal lines at 50 points equally spaced along 
the curve. We interpolated the intensity onto each of these lines, and located the peak 
emission. We used this sequence of points as an improved guess for the location of the spine, 
and iterated the procedure until no point along the spine moved by more than 0.1 pixel (1 
pixel = 3"). 

Our analysis "straightened" the filament by taking the radius to any point as the nearest 
perpendicular distance to the curve defining the spine. This procedure can be problematic 
for points on the concave side of a curved spine and more distant than its radius of curvature, 
which can, in principle, intersect with more than one orthogonal. However, Gil. 11-0. 12 is 
sufficiently linear that this was not a problem. Note that one of us used this procedure 
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Fig. 1. — The Gil. 11-0. 12 filamentary infrared dark cloud at 850 /mi. The labeled polygons 
are the regions chosen for analysis. Each of these regions contains one bright core (shown) 
that we removed to emphasize the filamentary structure. 
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previously in an analysis of the Orion integral-shaped filament (Matthews et al. 2001). 

Two regions selected for analysis are shown in Figure 1. The brightest peaks were 
clipped from the data as shown, since the contours within these regions are quite round 
and filament models do not apply. We binned the data to reduce the apparent scatter, and 
estimated the uncertainty from the standard deviation about the mean within each bin. We 
show scatter plots of the binned density as a function of radius for each region, along with 
fits to the data, in §6.1-6.3. 

In Paper I, we fit a non-pressure-truncated Ostriker model to each of two regions of 
Gil. 11 -0.12. We determined that the southern segment, corresponding roughly to our region 
1, could be fit by an Ostriker model with n c = (3.3 ±0.7) x 10 4 cm -3 and r = 0.10 ± 0.02 pc. 
We found that the northern segment, corresponding to our region 2 roughly, has a central 
number density of n c = (1.1 ± 0.3) x 10 4 cm -3 and a core radius of r = 0.12 ± 0.02 pc, 
assuming a dust temperature of 13 K and a dust emissivity of k$5o = 0.02 cm 2 g -1 . At 
the 3.6 kpc distance of Gil. 11-0. 12 the filament is only marginally resolved on these scales 
and thus care was taken to properly convolve the model to fit the data. The full width half 
maximum of an Ostriker filament is several r , however, and thus the core radius fits are 
reasonably well determined. 

4. Models and Observational Constraints 

Section 2.3 discussed how the poloidal and toroidal field components can compensate 
for each other to produce a 1-parameter set of models with approximately the same density 
(or surface density) structure. The degeneracy in T z and is approximate for the FP 
model, but the corresponding degeneracy in (3 and 9 is exact for the GS model. The goal 
of our analysis is to explore completely the parameter space of each model and to find the 
largest possible class of solutions that is consistent with the data. This determines the 
structure of the parameter space in detail, which would illuminate any other less obvious 
model degeneracies that might exist. 

It is useful to enumerate the parameters of each of the 3 models discussed in §2 be- 
fore discussing the results of our parameter space search. There is a single dimensionless 
theoretical parameter C for the truncated Ostriker model, and 3 dimensionless theoretical 
parameters for each of the magnetic models: C,f3,9 for the GS models and C, r 2 ,r^ for 
the FP model. All models require 5 additional parameters that reflect various observational 
complexities and unknown quantities. Two dimensional parameters p c and a determine the 
density scale and core radius. One parameter % encodes the unknown inclination of the 



filament relative to the plane of the sky. Finally, 2 additional dimensional parameters are 
required to convert the surface density £ calculated from the models to a predicted 850yum 
flux density / that can be compared with the observations. These are the flux conversion 
factor FCF, defined by equation 20 below, and a DC flux level foe that accounts for the 
unknown zero point of the 850/xm map. The latter parameter is necessary because the JCMT 
measures only a differential signal between chop positions and is insensitive therefore to the 
absolute zero point of the flux in the map. The total number of parameters is 6 for the 
Ostriker model and 8 for the magnetic models, but all 5 of the dimensional parameters are 
constrained by independent observations or reasonable assumptions. If one is concerned with 
the density structure only, and not with the magnetic field geometry, then the degenerate 
parameters (3 and 9 of the GS model can be replaced by a single independent parameter 
k (see discussion following equation 10), thereby decreasing the parameter space of the GS 
model to 7 dimensions. 

We determined the velocity dispersion within the filament from C 18 (2-1) and 13 CO (2- 
1) maps of the southern half of the Gil. 11-0. 12, including all of region 1. These observations 
were obtained with RxA3 on the JCMT during several observing sessions in 2000 July. 
Typical line widths were around 0.9 km s _1 for C 18 a 1.2 km s _1 for 13 CO. Since the 
C 18 data were often quite noisy, we loosened the constraint on the 1-dimensional velocity 
dispersion of the gas to be : 

0.8 km s" 1 < a < 1.2 km s _1 . (18) 

We assumed that the central density along the spine of the filament is in the range 

10- 21 g cm" 3 <p c < 10~ 17 g cm" 3 , (19) 

which corresponds to the range 2.6 x (10 2 — 10 6 ) cm- 3 in number density. The central density 
actually determined by our modeling procedure is well inside of these bounds. 

The flux conversion factor FCF is given in terms of the 850/im extinction coefficient 
«85o and dust temperature Td as follows: 

(20) 

(see Paper I). We combine the observational constraints on Td and k 85 o into a single constraint 
on FCF, since the temperature and extinction coefficient enter into our models only through 
FCF. The bounds on the Td and k 850 are given by 

8K < T d <15K 
0.01 cm 2 g- 1 < k 850 < 0.02 cm 2 g"\ (21) 
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which together imply the constraint 

0.27 g cm' 2 Jy' 1 beam < FCF < 1.92 g cm' 2 Jy' 1 beam. (22) 

The assumptions that lead to these constraints are discussed in Paper I. Note that the Os- 
triker solutions presented in Paper I required a k 85Q value of 0.02 cm 2 g' 1 and a temperature 
of 13 K. These values are at or near the upper bounds given by equation 21. 

The DC flux level foe was assumed to be small compared to the flux level over most of 
the map: 

\f DC \ < 0.25 Jy beam' 1 . (23) 

We verified that the actual fnc values obtained from the model fitting were well within this 
range. Thus, we are confident that we have not been too restrictive in setting this constraint. 

Finally, we assumed for all models that the axis of the filament lies with 45 degrees of 
the plane of the sky: 

-\<i< I m 

Note, however, that this quantity is unconstrained by observations, and it is conceivable that 
Gil. 11-0. 12 might be even more severely inclined. 



5. Parameter Space Exploration and Data Modeling 

5.1. Genetic Algorithm Basics 

Biological evolution can be thought of as a powerful and robust search and optimization 
algorithm. Evolution drives individuals within a population toward better adapted forms 
characterized by greater fitness and improved survivability. Genetic algorithms are a class 
of search and optimization algorithms that are based on biological evolution, in which the 
parameters of a model are encoded on the "genome" of an individual, and the goodness 
of the resulting solution translates into the fitness of the individual, and hence determines 
the likelihood that it will survive. Genetic algorithms are often associated with compu- 
tational intelligence research in computer science, since they learn by experimenting with 
their computational environment, and are capable of finding innovative solutions to complex 
optimization and design problems with little human intervention. 

Traditional genetic algorithms encode parameters on "genes" composed of low-cardinality 
(often binary) integer strings and process information via crossover, mutation and selection 
operators. The biological counterparts of these genetic operators motivate their respective 
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roles. Crossover combines the parameters represented by two "parent" individuals to pro- 
duce "offspring" with intermediate characteristics. Mutation provides occasional random 
perturbations to a small subset of the population. The selection operator serves as an infor- 
mation filter that destroys poor solutions, while allowing good solutions to survive, often by 
a tournament-based ranking scheme. These processes operate on a population of solutions 
iteratively, resulting in a global improvement in the quality of the solutions with each succes- 
sive generation. Goldberg (1989) is the standard introduction to the subject and provides 
code for a basic genetic algorithm. 

Genetic algorithms are guided stochastic search algorithms, which are much more ef- 
ficient than random search techniques. It is the non-random aspects of evolution, namely 
natural selection and crossover, which make evolution feasible and powerful as an algorithm 
for search and optimization. In general, genetic algorithms are extremely robust, and perform 
well in the presence of multiple local optima, discontinuities, noise and functions with de- 
ceptive gradients. They are suitable for mathematically difficult problems, high-dimensional 
parameter space searches, problems that aim to find multiple optima or entire classes of 
solutions, and multi-objective problems. However, they are not efficient for easy problems 
that can be handled by straightforward calculus-based techniques. 

One of us (J.F.) has developed a multi-objective genetic algorithm called "Ferret" with 
features that are optimized for parameter exploration and optimization problems encountered 
in the physical sciences. The algorithm has some similarities to the multi-objective algorithm 
discussed by Fonseca & Fleming (1993), but we have incorporated several new features to 
improve the flexibility and performance of the code. A detailed description of the algorithm 
is beyond the scope of this paper, but the most significant challenges involved in developing 
the code are discussed below. 



5.2. Solution Class-Finding & Multi-Objective Genetic Algorithms 

The data-fitting problem discussed in this paper is an example of a solution "class- 
finding" problem, where the goal is to find a class of acceptable solutions to a problem that 
obeys a user-defined acceptance criterion. For the Gil. 11-0. 12 problem, we define a x 2 that 
is representative of a good solution, and search for the class of all solutions with x 2 less than 
this value. 

Ferret is designed with multi-objective capabilities that will be useful for future polariza- 
tion modeling work, which would fit 850/iin continuum, the polarization position angle, and 
the polarization percentage simultaneously. It is important to note that most of the computa- 
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tional machinery needed for a multi-objective code is also required for solution class-finding. 
Ferret's multi-objective capabilities and the details of multi-objective optimization are dis- 
cussed in Fiege et al. (2004), where the code is being used model the internal structure of 
the Jovian moon Europa, using observational constraints obtained from the Galileo mission 
to Jupiter. 

Simple genetic algorithms suffer generally from a property known as "genetic drift", 
which causes solutions to converge toward a single arbitrary solution within the optimal set. 
Genetic drift is analogous to an undesirable loss of diversity within a biological population. 
This is problematic because it artificially reduces the extent of the solution class discovered 
by the algorithm. This problem must be treated carefully because it is very difficult for 
most genetic algorithms to recover, once diversity is lost in the solution class. Sophisticated 
countermeasures are required to combat genetic drift, which make class-finding and multi- 
objective codes much more complex than a simple genetic algorithm. Ferret uses a niching 
technique similar to that used by Fonseca & Fleming (1993), in addition to several other new 
techniques that are unique to our code. Ferret is the only multi-objective genetic algorithm 
that we are aware of that can recover completely from an extreme (artificially induced) loss 
of solution diversity within a few generations. 

Genetic algorithms have not been used widely in astrophysics, even though they are well- 
suited for many problems of astrophysical interest. Charbonneau (1995) provided a basic 
single-objective genetic algorithm PIKAIA and used it to solve several illustrative problems 
including modeling the rotation curves of galaxies, time series analysis of Doppler velocities 
in spectral lines of 5 Scuti stars, and modeling spherically symmetric magnetized winds 
from solar-type stars. This paper and the code provided therein provide a good alternate 
introduction to genetic algorithms. 

5.3. Code Testing and Verification 

We tested Ferret on a variety of multi-parameter multi-objective test problems of varying 
complexity. These tests included the following: 

1. Simple functions with multiple minima 

2. Noise tolerance tests 

3. Functions with deceptive gradients, which often fool calculus-based techniques 

4. Mult i- dimensional solution class-finding problems, using niching and other diversity- 
preserving strategies 
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5. Class-finding problems with discontinuities and complex geometries 

6. Multi-objective problems, using niching and other diversity-preserving strategies 

7. Discrete and semi-discrete problems where some of the variables take only discrete 
values 

8. Critical-parameter detection problems, where the solutions reside in a low-dimensional 
subspace of the parameter space. The code improves performance by determining 
the relative importance of each parameter and focusing on the parameters that are 
significant to the problem. This allows solutions to be found in test problems designed 
with up to 500 parameters. 

The code performed well and located the solution set for each of these problems. 

We show the results for a test problem that is similar in spirit to the parameter search 
conducted in this paper. We consider an N-dimensional function f — fi f 2 , where f\ and f 2 
are defined as follows: 



fi = sign 

N 



2nr \ f 2nr 

cos — — — cos 6 + 



Ai J V Ai 



f 2 = g 2 r und(l^^)-/Z«^(!^5l) 



(25) 



where r are 9 are cyclindrical coordinates in the (xi,x 2 ) subspace defined by 

r = sj (x\ + xj) 
6 = tan' 1 (x 2 /x 1 ), (26) 

and r is a constant. The functions "sign(x)", "round(x)" and "floor(x)" are the familiar 
discontinuous functions that are available in most programming languages. We designed 
this rather complicated function as a non-trivial test of several features of the code that are 
critical to its success on data modeling problems. Specifically, we test the ability of the code 
to work efficiently in very high-dimensional parameter spaces, to avoid being trapped by 
numerous local minima, and to function as a solution class-finder by distributing solutions 
over a finite-volumed optimal region. Furthermore, our test problem has a simple solution, 
but cannot be solved by any calculus-based numerical approach because the gradient is zero 
everywhere, except at discontinuities. 

The function f\ is a discontinuous function with a long, flat-bottomed spiral whose 
minimum has thickness 2*r and repeats with radial wavelength Ai in the (xi,x 2 ) subspace. 



-16- 



In this finite domain, fi breaks up into 9 distinct minima of equal depth. The function 
f 2 is another discontinuous function designed to generate many local minima, spaced by 
wavelength A2, and a single global minimum of width A2 centred at Xq. We chose N = 20, 
Xq = 5, and a range of —10 < Xi < 10 for all parameters. Figure 2a shows one term of f 2 , 
where A2 = 1 and, over the same range as above. Figure 2b shows the sum over any two 
terms, where we have zoomed in to the region near the global minimum at (5,5). 

On our chosen range, each term of f 2 generates 19 local minima, including a single 
global minimum. The summation over N — 2 terms results in 19 N ~ 2 minima. Multiplying 
by the 9 minima in f\, we find that the combined function / has 9 x 19 N ~ 2 = 9.4 x 10 23 
isolated, finite-volumed minima, of which only 9 are global minima. Our goal is to map out 
the 9 global minima while evading all other local minima, which should reveal the spiral 
structure of f\ in the (xi,x 2 ) projection, and look like a 1 x 1 square box centred on (5,5) 
in any projection that does not include x\ or x 2 . 

We ran the code using a population size of 1000 for 500 generations, and performed 
several trials. In each trial, all global minima were found and covered with points in less 
than 100 generations. Subsequent generations increased the density of points and improved 
the resolution of the boundaries. Figure 2c shows the projection of the solution set in the 
(xi, x 2 ) subspace, and Figure 2d shows an example of another projection. All subspaces that 
do not include x± or x 2 look like Figure 2d. 

This problem is more difficult that the parameter searches required to model the Gil. 11- 
0.12 cloud because our astrophysical models have fewer (6-8) parameters, and few local 
minima, other than the global minimum. We conclude that our code succeeded on this test 
problem, and should be adequate for the modeling conducted later in this paper. 

5.4. The Objective Function for Gil. 11-0. 12 models 

In Paper I, we fit two Ostriker (1964) solutions to the 850/im continuum data for Gil. 11- 
0.12. We extend this analysis by exploring and constraining the parameter spaces for the 3 
models discussed in §2. To accomplish this goal, we need to locate the class of solutions for 
each model, which has a reduced x 2 < Xmin where Xmin > 1 but s ^ m °f order unity. We also 
need to rule out solutions or particular regimes of each model (i.e. poloidal-field dominated 
models, toroidal-field dominated models, etc.) where this condition is not met. 

We verified the success of the code by checking the values of x 2 a posteriori. We ran 
the genetic algorithm typically for 500 generations, although a stable population of models 
satisfying our x 2 criterion usually developed after only ~ 30 — 40 generations. Subsequent 
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Fig. 2. — High dimensional test problem for Ferret genetic algorithm with many (~ 9.4 x 10 23 ) 
local minima. Panel a (top left): One term in the function f 2 given by equation 25. Panel b 
(top right): Two terms of function f 2 , zooming in on the global minimum at (5,5). Panel c 
(lower left): Projection of the solution set in the (x±, x 2 ) subspace showing the flat-bottomed 
spiral minimum given by function f±. Panel d (lower right): Projection of the solution set 
in the (x 7 ,x 20 ) subspace. The solution set looks like this figure in all subspaces that do not 
include X\ or x 2 . 
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generations served to fill out the optimal class of solutions and to improve the resolution of 
the edges of the optimal region of parameter space. In practice, we save all solutions that 
are in the optimal set of each generation, and filter through all generations at the end of the 
run to select all models from all previous generations that satisfied the x 2 criterion. Our 
results show that the best solutions discovered are usually surrounded by a sheath of poorer 
solutions, as one might expect. We ran the code several times while varying the population 
size and other genetic algorithm control parameters, to ensure that that the solution set was 
repeatable. 

We normally display the output from our genetic algorithm as a scatter plot that repre- 
sents a 2-dimensional projection of the region of parameter space that agrees with the data, 
according to our x 2 criterion. In practice, a well-designed multi-objective or class-finding 
genetic algorithm distributes solutions as evenly and completely as possible throughout the 
optimal region. However, there is no strict mathematical guarantee that this will occur. The 
safest interpretation of our scatter plots is simply that the locus of points covers the multi- 
dimensional region of parameter space allowed by the observations, without attempting any 
further inference regarding the uniformity of the points or the probability distribution of the 
solutions. 

As noted in §5.2, multi-objective problems require nearly the same computational ma- 
chinery as the single-objective class-finding problem that is the focus of this paper. Future 
work will focus on 2-dimensional polarization modeling, which will fit the intensity map, the 
polarization position angle and the polarization percentage simultaneously across the map. 

6. Results 

6.1. The Truncated Ostriker Model 

Figure 3 shows examples of truncated Ostriker model fits to regions 1 (left panel) and 
2 (right panel). The radial profile of region 1 is smoother and less noisy than that of 
region 2, which results in better constrained fits. The grayscale indicates the quality of the 
solutions: black designates solutions with x 2 < 1.167, medium gray designates solutions with 
1.167 < x 2 < 1-333, and light gray is for solutions with 1.333 < x 2 < 1.5. We discuss the 
distribution of these solutions within the parameter space of the Ostriker model below. 
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Fig. 3. — Sample Ostriker model fits to the data. Region 1 is on the left and region 2 
is on the right. The grayscale indicates the quality of the fit: black indicates solutions 
with x 2 < 1-167, medium gray indicates solutions with 1.167 < x 2 < 1.333, and light gray 
indicates solutions with 1.333 < x 2 < 1.5. We use this convention in all of our scatter plots, 
except for those of §7. We always plot points in order of decreasing range of x 2 , so that 
darker points, representing solutions of higher quality, are plotted after lighter points and 
may obscure them in some cases. 
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6.1.1. Observationally Constrained Parameters 

Figure 4 shows the set of truncated Ostriker solutions in the observationally constrained 
part of the parameter space (p c , a, FCF) for region 1. We found solutions with central 
number densities in the range 0.98 x 10 4 cm -3 < n c < 2.4 x 10 4 cm -3 and velocity disper- 
sion constrained to the range 0.86 km s" 1 < a < 1.2 km s" 1 . The flux conversion factor 
is constrained to values near the lower limit of the allowed range given by equation 22: 
0.27 g cm" 2 Jy" 1 beam < FCF < 0.53 g cm" 2 Jy" 1 beam. 

Our analysis found a central number density that is somewhat lower than the value 
given for the southern region by Paper I, which gave n c = (3.4 ±0.6) x 10 4 cm -3 for region I. 
The grayscale indicates the quality of the solutions, following the same convention outlined 
in the caption of Figure 3. There are no apparent gaps in the solution set, and the best 
solutions are surrounded by a sheath of poorer solutions, as expected. 

We obtained similar results for region 2, but the solutions are much less well-constrained 
than in region 1. The central number density of region 2 is constrained to the range 3.5 x 
10 3 cm~ 3 < n c < 2.0 x 10 4 cm~ 3 . The velocity dispersion spanned the entire range of 
allowed values from 0.8 km s~ x < a < 1.2 km s -1 . Allowed values of the flux conversion 
factor ranged from 0.31 g cm~ 2 Jy _1 beam < FCF < 1.6 g cm~ 2 Jy~ x beam, which covers 
a significant fraction of the allowed range. The central number density is consistent with the 
value of (1.1 ± 0.3) x 10 4 cm -3 obtained in Paper I, but the range of densities that we found 
indicates that the uncertainty in their estimate was probably underestimated. 

As we discussed in §4, our modeling technique combines the temperature (T) and 850 
/iin emissivity (/«85o) into a single flux conversion factor (FCF). Each value of the FCF 
corresponds to a one parameter set of (T, /tsso) values, according to equation 20. Figure 5 
plots a random sample of these curves to show the possible ranges of these parameters for 
region 1 (left panel) and region 2 (right panel). The dashed lines in this figure indicate the 
observational constraints on the temperature. The results for region 1 require values of T 
and k 850 that are in the upper range of the values given by equation 21. This is consistent 
with the model presented in Paper I, which used a temperature of 13 K and an emissivity 
of 

K 850 — 0.02 cm g . Region 2 is less well-constrained, but it is consistent with a broad 
range of T and k 850 covering most of the allowed range. 

The solution sets for both regions span the entire range of inclination angles allowed 
by equation 24. We did not find any interesting correlations between the inclination angle 
and any other quantity. Obviously, the inclination angle must be correlated with other 
parameters, but the underlying relationship is not evident in any 2-dimensional projection 
of the parameter space that can be displayed here. 



- 21 - 




n (10 4 crrf 3 ) n (10 4 cm' 3 ) 

C C 

Fig. 4. — Two views of the observationally constrained part of the parameter space 
(p c ,a,FCF) of the truncated Ostriker models for region 1, with \ 2 < 1-5. The grayscale 
indicates the quality of the solutions, according to the convention of Figure 3. The best 
solutions form a core surrounded by an envelope of poorer solutions. 




Fig. 5. — Constraints on the temperature and 850/im extinction coefficient of the truncated 
Ostriker model. We show a random sample of solutions for region 1 (left) and region 2 
(right). Each curve corresponds to a single value of the FCF for a single model, as defined 
by equation 20. The grayscale follows the convention of Figure 3. 
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6.1.2. Theoretical Parameters 

Figure 6 shows the virial parameters of the truncated Ostriker model solutions. The 
quantities m/m V i r and Ps/ < P > obey equation 4 for the truncated Ostriker model, and 
equation 15 for more general magnetic models. FP1 showed that these virial quantities 
provide an important constraint on the geometry of the magnetic field, since equation 15 
allows one to deduce the value of the magnetic parameter .M/|W| if m/m vir and Ps/ < 
P > can be estimated. Truncated Ostriker models always lie on the magnetic neutral line 
(dotted line) of this diagram. It is notable that Gil. 11-0.12 appears to have a lower value 
of Ps/ < P > and a much higher value of m/m vir than the filaments that were studied by 
FP1. The fact that good fits were found to the Ostriker model indicates that the line mass of 
Gil. 11-0. 12 is high enough for self-gravity to hold the filament together without additional 
binding energy due to a toroidally dominated magnetic field. 

6.2. The Generalized Stodolkiewicz Model 

Figure 7 shows examples of GS model fits to our data. Both of the regions that were 
analyzed are well-fit by the model, but the solutions are not as well-constrained in region 2. 

6.2.1. Observationally Constrained Parameters 

Figure 8 shows the distribution of GS models for region 1, in the observationally con- 
strained part of the parameter space given by (p c ,a, FCF). The set of models that fit the 
data are constrained approximately to a curved 2-dimensional sheet in this subspace. Good 
fits are available for the entire range of a and FCF explored. The GS model is consistent 
with a range of central number density such that 9.8 x 10 3 cm -3 < n c < 1.3 x 10 5 cm -3 . 

Models were less well-constrained for region 2, and the corresponding figures are omitted 
for brevity. Region 2 is consistent with central number densities of 4.0 x 10 3 cm -3 < n c < 
3.5 x 10 4 cm -3 , which are a factor of 2.5 — 3.5 lower than region 1. The central number 
density range is consistent with the values that we obtained for the truncated Ostriker model, 
but the GS model has a higher upper limit. 

The solution set for the Stodolkiewicz model spans the entire range of inclination angles 
allowed by equation 24. There are no obvious correlations between the inclination angle and 
any other quantity, for reasons that are discussed at the bottom of §6.1.1. 
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Fig. 6. — The virial parameters of the truncated Ostriker model solutions. Results for region 
1 are shown on the left, and region 2 on the right. The quantities m/m vir and Pg/ < P > 
obey equation 4 for the truncated Ostriker model, and equation 15 for more general magnetic 
models. Magnetic models dominated by a toroidal field would lie to the left of this curve, 
and poloidal-field dominated models would lie to the right. The numbers labeling the curves 
indicate values of the magnetic parameter .M/|W| ( see equation 15). The gray box indicates 
the location of the toroidal-field dominated filamentary clouds studied by FP1. This box 
is almost entirely on the A^/|W| < side of the diagram. This region is inaccessible to 
the non-magnetic truncated Ostriker solutions, which are always constrained to lie on the 
magnetically neutral line (heavy line). The grayscale follows the convention of Figure 3. 
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Fig. 7. — Sample GS model fits to the 850 /im intensity data. Region 1 is on the left and 
region 2 is on the right. The grayscale indicates the goodness of fit, following the convention 
of Figure 3. 




Fig. 8. — Two views of the observationally constrained part of the parameter space 
(p c ,cr,FCF) of the GS model for region 1, with x 2 < 1-5- The allowed models reside 
approximately within a curved 2-dimensional sheet in this subspace. The grayscale indicates 
the quality of the solutions, according to the convention of Figure 3. 
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6.2.2. Theoretical Parameters 

The theoretical parameters for the GS model are the concentration parameter C (equa- 
tion 6) and the magnetic parameter k, discussed earlier in this section. There were no 
interesting correlations between these quantities. For region 1, solutions filled the rectangu- 
lar region of parameter space decribed by 0.89 < C < 2.0 and 0.1 < k < 12. For region 2, 
solutions filled the rectangular region of parameter space decribed by 0.75 < C < 2.0 and 
0.1 < k < 5.1. 

Figure 9 shows the virial diagrams for the allowed GS solutions that are consistent with 
each of the two regions that were studied. Both regions are consistent with a poloidally 
dominant magnetic field with 7W/|W| < 1. Note, however, that toroidal-field dominated 
solutions are not possible within the context of the GS model. 

Figure 10 shows how the magnetic energy ratio M/K scales with the magnetic parameter 
(1 + cos 2 6)/2{3 within the allowed range of solutions for region 1. The plasma (3 parameter 
is equal to the inverse of the magnetic energy ratio, and the scatter in this figure is due to 
the range of magnetic pitch angles 9 allowed by the GS model. 

6.3. The Fiege & Pudritz Model 

Figure 11 shows examples of FP model fits to both regions of Gil. 11-0. 12 that were 
analyzed. Both regions are well-fit by the model, but the solutions are less well-constrained 
in region 2. This section discusses the distribution of these solutions in the parameter space 
of the FP model. 



6.3.1. Observationally Constrained Parameters 

Figure 12 shows two different viewpoints of the distribution of FP models for region 1, 
in the observationally constrained part of the parameter space given by (p c ,a,FCF). The 
set of models that fit the data occupy an approximately planar region in this subspace. We 
obtained good fits for essentially the entire range of a and FCF explored. The FP model is 
consistent with central number densities in the range 9.1 x 10 3 cm~ 3 < n r < 1.1 x 10 5 cm -3 . 
The highest number density that we obtained for an Ostriker model was 2.4 x 10 4 cm -3 ; 
thus, the FP model is consistent with central densities of up to ~ 5 times higher. This 
is significant because higher densities may be more consistent with the critical densities of 
molecules such as H 2 CO (Carey et al. 1998), HCO+ and CS (Redman et al. 2000), and 
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Fig. 9. — The virial diagram for the GS model. Region 1 is shown on the left and region 2 
on the right. The diagram is explained by the caption of Figure 6, and the grayscale follows 
the convention of Figure 3. Note that the magnetic field of the GS model always contributes 
to net magnetic support. Thus the A^/|W| < (magnetic confinement) side of the diagram, 
where the filaments studied by FP mostly reside, is inaccessible to this model. 




Fig. 10. — The magnetic energy ratio M/K of the GS model as a function of the magnetic 
parameter k = (1 + cos 2 9)/2(3. Only solutions that are consistent with region 1 are shown. 
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N 2 H + (Feldman et al. 2003). Higher densities are also preferred by current ideas on dust 
coagulation, which assume densities of ~ 10 6 cm -3 , as discussed in Paper I. The velocity 
dispersion and FCF spanned the full range allowed by the search. 

We found similar results for region 2 within the following parameter ranges: 2.7 x 
10 3 cm" 3 < n c < 4.0 x 10 4 cm -3 . The number density range for both regions is similar to 
our results for the GS model, but much higher than values obtained for the Ostriker model. 

Figure 13 shows the allowed values of T and k 850 for a sample of FP models for region 1 
(left) and region 2 (right). FP solutions exist for nearly any value within the range allowed 
by equation 21. 

The solution sets for both regions span the entire range of inclination angles allowed by 
equation 24. We did not find any obvious correlations between the inclination angle and any 
other quantity for the FP model. 

6.3.2. Theoretical Parameters 

Figure 14 shows the distribution of FP model solutions in the theoretical part of the 
parameter space (r z , T^, C) for region I. This region is consistent with values of T z < 15 
and < 3.7. There is a very clear trade-off between flux to mass ratios T z and T^, which 
lie approximately on a curve in the (T z , T^) subspace, such that increases with increasing 
T z . This trend occurs because T z and have an opposite effect on the radial structure of 
the cloud. Equation 17 and the discussion of §2.3 show that any value of the magnetic virial 
parameter At/|W| corresponds to a one parameter set of T z and T^. Figure 14 demonstrates 
that this exact degeneracy in the virial balance is also manifest as an approximate degeneracy 
in the radial density structure of the cloud. This degeneracy is expected from the results 
obtained by FP1, who showed that poloidal and toroidal fields have opposite effects on the 
shape of the density profile. Dominant toroidal fields result in centrally peaked density 
profiles with a slow decrease outside of the core region, while dominant poloidal fields result 
in flat central density profiles with a sharp decrease in the outer regions. Our results are an 
indication that these two effects can be critically balanced against each other to explain a 
given density profile for a significant range of Y z and T^. 

It is very interesting that the swath of solutions shown in Figure 14 just grazes the 
non-magnetic limit (r z = 0, = 0), which is equivalent to a truncated Ostriker model. 
The Ostriker model is consistent with the data, but it can be viewed as an outlying model 
near the boundary of a more general class of magnetic solutions. 



-28- 




r (pc) r (pc) 



Fig. 11. — Sample Fiege & Pudritz model fits to the 850 /zm intensity data. Region 1 is on 
the left and region 2 is on the right. The grayscale indicates the goodness of fit, following 
the convention of Figure 3. 




Fig. 12. — Two views of the observationally constrained part of the parameter space 
(p c ,a, FCF) of the FP model for region 1, with \ 2 < 1-5- The grayscale indicates the 
quality of the solutions, according to the convention of Figure 3. 
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Region 1 is consistent with values of the concentration parameter C in the range 0.89 < 
C < 2, which corresponds to a truncation range of 7.8 < Rs/tq < 100. We found that Rs is 
poorly determined within this range. This is not surprising because the density is very low 
at these large radii and has little effect on the surface density structure. The concentration 
parameter is slightly correlated with T z and tends to slowly decrease as T z increases. This is 
explained by the opposite roles of these parameters in the surface density structure. Solutions 
with large values of Y z have large, flat central density profiles with more emitting material 
near the center than low T z solutions. Decreasing C compensates for high T z by decreasing 
the column depth of the surrounding gas. 

The corresponding figures are shown in Figure 15 for region 2. Region 2 is consistent 
with values of T z < 15 and < 6. and T z are correlated, as in region 1, but less tightly. 
The solutions have a similar range of concentrations as in region 1: 0.80 < C < 2.0. 

FP defined a magnetic energy parameter M/K to determine the relative importance 
of the volume-averaged magnetic energy density M, compared with the volume-averaged 
turbulent energy density K. These quantities are defined as follows: 

= ! v {Bl + Bl)dV 
8irV 

3 2 

= 2 <p><a >J 

where all quantities are defined as in §2.3. FP considered M/K pa 5 to be the upper limit 
for reasonable models, and 0.5 < M/K < 2 to be the preferred range. Figure 16 shows the 
distribution of M/K for region 1. This figure shows that M/K is almost exactly a function 
of Y z alone, and that solutions with T z < 15 have M/K < 4. Figure 17 is the corresponding 
one for region 2. There is more scatter in this figure, but M/K remains highly correlated 
with T z . Solutions with T z < 15 are consistent with M/K < 5. 

Figure 18 shows the virial parameters of the FP solutions for regions 1 (left panel) and 
2 (right panel). We find that .M/|W| > for all models or Region 1 and most models 
of Region 2, which indicates that most models are consistent with a dominant poloidal 
magnetic field. Some magnetically neutral models A^/|W| = are allowed 2 , and the few 
toroidal field-dominated models that are allowed are characterized by rather weak fields. 
This is strikingly different from the results obtained by FP1, which are indicated by the 
shaded box. Toroidal-field dominated models, which were the focus of FP1, are unlikely 



2 Magnetically neutral models are not necessarily non-magnetic. Filaments with M/|W| = can be 
obtained by critically balancing the toroidal field and poloidal field terms in equation 17. 
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Fig. 13. — Constraints on the temperature and 850/iin extinction coefficient of the FP model. 
We show a random sample of solutions for region 1 (left) and region 2 (right). Each curve 
corresponds to a single value of the FCF for a single model, as defined by equation 20. The 
grayscale follows the convention of Figure 3. 




Fig. 14. — Distribution of the theoretical parameters of the FP model (T z , T^, C) for region 
1. (left) T z and are tightly correlated and lie approximately on a curve in the (T z , 
T^) subspace. (right) The lower bound of the concentration parameter C tends to increase 
slightly as T z increases. 
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Fig. 15. — Distribution of the theoretical parameters of the FP model (T z , T^, C) for region 
2. (left) T z and are correlated and lie approximately on a curve in the (T z , T^) subspace. 
(right) The concentration parameter C. 
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Fig. 16. — The FP magnetic energy parameter for region 1. (panel a., left) We show the 
correlation of M/K, defined by equation 28, with T z . (panel b., right) We overlay several 
values of M/K on a graph of vs. T z . The dashed lines represent the standard deviation 
of the value of T z corresponding to each value of M/K. 
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for this filament, and our results indicate that Gil. 11-0. 12 is an excellent candidate for a 
magnetically supported filament. 

7. Polarization Predictions 

Fiege & Pudritz (2000c) developed a technique to simulate submillimetre polarization 
maps of filamentary cloud models. They considered only the FP model, but we generalize 
their analysis to predict the polarization structure due to all GS and FP solutions that 
are consistent with the radial structure of the Gil. 11-0. 12 IRDC. The Fiege & Pudritz 
(2000c) formalism combines the uncertainties due to poorly contrained dust properties and 
grain alignment mechanisms into a single parameter < a >. This parameter determines 
the maximum polarization percentage that can be obtained from a region threaded by a 
magnetic field with the optimal geometry, which corresponds to a uniform magnetic field 
parallel to the plane of the sky. We assume that < a >= 0.15, which leads to a maximum 
polarization percentage of 15.4% by equation 17 in Fiege & Pudritz (2000c). 

Fiege & Pudritz (2000c) found that filaments aligned parallel to the plane of the sky 
are consistent with three general types of polarization patterns. Type 1 patterns are charac- 
terized by polarization vectors that are aligned exactly parallel to the axis of the filament, 
and are consistent with toroidal-field dominated models. Type 2 patterns have polarization 
vectors that are aligned exactly perpendicular to the filament, and are due to poloidal- field 
dominated models. Type 3 patterns occur in models where the poloidal and toroidal fields 
are of comparable strength. Type 3 polarization maps are aligned perpendicular to the 
filament in some regions (especially near the poloidally dominated core), but are aligned 
parallel to the filament in other regions (i.e. a toroidally dominated envelope). Matthews 
et al. (2001) used this technique to model the 'integral-shaped' filament in Orion A, which 
they found to be consistent with a Type 1 polarization pattern. Fiege & Pudritz (2000c) 
also showed that Type 3 patterns are most likely for FP filaments that are inclined relative 
to the plane of the sky. 

We used this scheme to classify the polarization patterns of all GS and FP models 
computed in this paper. All GS solutions that are consistent with the data are Type 2 
or Type 3 polarization patterns. A polarization map that shows a Type 1 pattern would 
rule out the GS model. Figure 19 shows that the maximum polarization percentage may 
range from about 4% to 15% for both regions. This figure represents the three x 2 bins used 
elsewhere in the paper, plus the 3 possible polarization types, of which only 2 actually occur. 

The left panel of Figure 20 shows the distribution of the type of pattern with respect to 
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Fig. 18. — Virial diagram for FP solutions for regions 1 (left) and 2 (right). All quantities 
are as defined in Figure 6. The FP model is the only model of the three models studied 
that permits solutions with A^/|W| < (net magnetic confinement). However, all solutions 
for region 1 are magnetically neutral or dominated by the poloidal field. Most solutions for 
region 2 are also dominated by the poloidal field. However, a few toroidal-field dominated 
solutions are permitted when the Ps/ < P > is in the upper range shown in the figure. 
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Fig. 19. — Submillimetre polarization predictions for the GS model. Region 1 is represented 
on the left, and Region 2 on the right. We classify the polarization structure according 
to the classification scheme of Fiege & Pudritz (2000c) and show the maximum predicted 
polarization percentage. All solutions are Type 2 (poloidal-field dominated) and Type 3 
(mixed) patterns. The colour scheme is consistent with Figures 20 and 21 below, which are 
the corresponding figures for the FP model. Points were layered as follows. Successively 
better (lower) x 2 bins always overlay poorer bins, as in previous figures. Within each bin, 
points were lain down in the order Type 1, followed by Type 2, and finally Type 3. Some 
overshadowing of points was inevitible between layers. 
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T z and for the FP model, applied to region 1. We find that Type 3 (mixed type) polariza- 
tion patterns occupy the largest part of the allowed parameter space. Type 2 (poloidal-field 
dominated) patterns occupy a smaller fraction of the parameter space, located at moderate 
T z and low T^. It is interesting that we also found a few Type 1 patterns, which appear 
to be toroidally dominated. These solutions are slightly dominated by their poloidal fields 
in a virial sense, or magnetically neutral, but would appear to be toroidally dominated in 
a polarization map. This is possible because the integrals that enter the virial analysis are 
quite different from the Stokes integrals. All of these solutions have very low values of T z 
and r^,. 

The right panel of Figure 20 shows how the maximum polarization percentage varies 
with the ratio r^/r 2 . Maximum polarization percentages ranging from 2-15% are consistent 
with the FP model for this region. A polarization map of Gil. 11-0. 12 that shows a Type 3 
pattern with maximum polarization percentage below ~5%, or any Type 1 pattern, would 
be consistent with the FP model, but not the GS model. It would be difficult to distinguish 
between the GS model and the FP model if a Type 2 pattern were observed, unless the 
peak polarization percentage was observed to be in the 2-4% range, since such low values 
are consistent with the FP model but not the GS model. However, detailed fitting to the 
actual polarization map might reveal other differences that help to distinguish between these 
models. 

It is clear that the type of polarization pattern, combined with the observed maximum 
polarization percentage provides an estimate of T^/T z . This is significant because r^/r 2 is 
the most important factor that determines the fragmentation timescale, the spacing of cores 
in the filament , and the masses of cores (Fiege & Pudritz 2000b). Therefore, polarization ob- 
servations of Gil. 11-0. 12, in combination with this work, would constrain the fragmentation 
properties of the cloud strongly, and determine the initial magnetohydrodynamic conditions 
that lead to star formation. 

Figure 21 is the corresponding figure for region 2. All three types of polarization patterns 
are possible for region 2. The polarization percentages have a similar range of allowed values 
as for Region 1, although the scatter is much greater. 

8. Discussion 

This paper is the first astronomical research problem attempted using a new multi- 
objective genetic algorithm "Ferret" designed by one of the authors (J.F.). We have tested 
this code thoroughly, and it performs extremely well on data-fitting problems where the goal 
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Fig. 20. — Submillimetre polarization predictions for FP solutions for region 1. We classify 
the polarization structure according to the classification scheme of Fiege & Pudritz (2000c): 
black dots indicate Type 1 patterns, medium gray dots designate Type 3 patterns and light 
gray is for Type 2 patterns. These patterns are explained in the text, and the layering 
of points is explained in Figure 19. [Note that this figure does not follow the grayscale 
convention used elsewhere in this paper, where the shade indicates the quality of the solution.] 
(panel a., left) The distribution of the polarization type with respect to T z and T^. (panel 
b., right) Polarization percentage versus r^/T 2 . 




Fig. 21. — Submillimetre polarization prediction of FP solutions for region 2. See Figure 20 
for explanation. 
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is to find the class of allowed solutions, rather than a single representative solution. This 
allows the user to understand the detailed structure of the parameter space, and gives an 
honest representation of the degeneracies of the models being tested. 

Our previous analysis of the Gil. 11-0. 12 infrared-dark cloud (Paper I) demonstrated 
that its density profile is steeper than the r~ 2 profile inferred from observations of other 
filamentary molecular clouds (Alves et al. 1998, Lada et al 1999, Johnstone & Bally 1999), 
and that the profile can be fit reasonably well by the non-magnetic Ostriker model. This 
paper provides a more detailed analysis that compares the radial structure of the Gil. 11- 
0.12 infrared-dark cloud with three models of filamentary clouds: the non-magnetic Ostriker 
model, the magnetic GS model, and the magnetic FP model. We found that the density 
profile of Gil. 11-0. 12 is consistent with the Ostriker model, but also with regimes of both 
magnetic models. Our technique allowed us to map out the allowed regions of parameter 
space thoroughly for each model. 

Our most significant result is that nearly all magnetic solutions that are in agreement 
with the data are dominated by the poloidal component of the magnetic field, or are mag- 
netically neutral. A few toroidal-field dominated solutions are allowed, but only for Region 
2, which has a much "noisier" radial structure than Region 1, and even these solutions are 
close to magnetic neutrality. Gil. 11-0. 12 appears to be very different from the filaments 
studied by FP1, which were inferred to have a shallow r~ 2 density profile and a dominant 
toroidal field component. Note that the FP model is the only model of the three models 
considered that can account for both toroidal-field dominated filaments and the poloidal-field 
dominated Gil. 11-0. 12 filament. 

We predicted the polarization patterns allowed by the GS and FP models of Gil. 11-0. 12 
that are allowed by the data. For the GS model, all of the polarization patterns are Type 2 
(poloidal-field dominated) or Type 3 (mixed) patterns Fiege & Pudritz (2000c). Thus, the 
GS model could be ruled out for the Gil. 11-0. 12 filament if future polarimetric observations 
reveal that this region is more consistent with a Type 1 (toroidal-field dominated) pattern. 

For the FP model, Type 3 patterns Fiege & Pudritz (2000c), which are consistent with a 
dominant poloidal magnetic field, occupy the greatest fraction of the parameter space. This 
is due primarily to the large range of inclination angles (-45° - 45°) present in the solution set, 
since Fiege & Pudritz (2000c) showed that type 3 patterns are the most common pattern 
when the filament is inclined relative to the plane of the sky. Type 2 patterns occupy a 
fraction of the parameter space, and occur at moderate values of T z and low values of T^. 
We found a few Type 1 patterns that would appear to be dominated by a toroidal field. Most 
of these solutions have very weak magnetic fields that are close to magnetic neutrality or 
slightly dominated by the poloidal field, according to the virial equation 15. We demonstrated 
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that the polarization class, combined with the maximum polarization percentage, provides an 
estimate of the ratio of magnetic flux to mass values Yz/Y^. This is significant because Fiege 
& Pudritz (2000b) demonstrated that this quantity is the dominant factor that controls 
the fragmentation timescale of a filament, the fragmentation wavelength, and ultimately the 
masses of clumps that form along the filament. Thus, a polarization map of Gil. 11-0. 12 
might constrain the magnetic field sufficiently well to provide an indication of the magnetic 
and dynamical conditions within the filament. Improved constraints could be obtained by 
detailed multi-objective modeling of the filament, which would fit the intensity structure, 
polarization vector orientations and polarization percentages simultaneously. 

Our polarization modeling also determined that the maximum polarization percentage 
may vary from 2-15% for the FP model, but is restricted to 4-15% for the GS model. Thus, 
observations of a Type 2 pattern with a low peak polarization percentage between 2-4% 
would not be consistent with the GS model. 

We have not considered non-axisymmetric magnetic field geometries in our analysis. It 
is interesting to note that a filament threaded by a straight magnetic field perpendicular to 
the filament axis would generate a Type 1 polarization pattern that may be similar to that 
of a toroidal field-dominated filament, if either the axis of the filament or the magnetic field 
is parallel to the plane of the sky. However, this is not true for the general case where neither 
the filament axis nor the field is parallel to the sky plane, because of projection effects. In 
this case, the projection of the field in the sky plane, and hence the polarization vectors, 
would have a contstant angle relative to the projected axis of the filament. The symmetry 
of a helical field results in polarization vectors that are always parallel or perpendicular to 
the projected filament axis, unless the filament is bent, in which case the polarization angle 
changes along the filament (Matthews et al. 2001; Matthews et al. 2002). A helical field 
cannot produce vectors that maintain a constant angle relative to the filament axis. Thus, 
it should be possible to distinguish between toroidal field-dominated helical field models and 
models with a non-axisymmetric orthogonal field for an ensemble of objects, although it 
may be difficult to distinguish between these possibilities for a single filament. It would be 
interesting to include non-axisymmetric equilibrium models in future analyses. However, no 
such models exist at present. 

Difficult data-modeling problems are encountered frequently in astronomy and physics, 
which may feature large parameter spaces, noisy data, multiple constraints, and multiple 
modeling objectives applied to heterogenous data sets. Our multi-objective genetic algorithm 
solves such problems reliably and efficiently. We hope that our code will prove useful for a 
variety of problems in the physical sciences. Our immediate goal is to build an automated 
system for multi-objective 2-dimensional modeling (i.e., 2 spatial map coordinates, rather 
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than a single radial distance) of continuum and polarization maps, which would allow us to 
model simultaneously the intensity, polarization position angle and polarization percentage 
for filaments and cores. In fact, we have already designed and tested a preliminary system 
that models the intensity and polarization data for one of the cores in the Orion Molecular 
Cloud J-shaped filament, by varying the magnetic geometry, the density structure, the 
temperature and dust emissivity, the position angle and inclination of the core relative to 
the plane of the sky This analysis will be the subject of future work. 
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